Chaotic Vibrations of the One-dimensional Wave Equation Due to a Self-excitation Boundary Condition Part I: Controlled Hysteresis

The study of nonlinear vibrations/oscillations in mechanical and electronic systems has always been an important research area. While important progress in the development of mathematical chaos theory has been made for finite dimensional second order nonlinear ODEs arising from nonlinear springs and electronic circuits, the state of understanding of chaotic vibrations for analogous infinite dimensional systems is still very incomplete. The 1-dimensional vibrating string satisfying wtt − wxx = 0 on the unit interval x ∈ (0, 1) is an infinite dimensional harmonic oscillator. Consider the boundary conditions: at the left end x = 0, the string is fixed, while at the right end x = 1, a nonlinear boundary condition wx = αwt − βw3 t , α, β > 0, takes effect. This nonlinear boundary condition behaves like a van der Pol oscillator, causing the total energy to rise and fall within certain bounds regularly or irregularly. We formulate the problem into an equivalent first order hyperbolic system, and use the method of characteristics to derive a nonlinear reflection relation caused by the nonlinear boundary condition. Since the solution of the first order hyperbolic system depends completely on this nonlinear relation and its iterates, the problem is reduced to a discrete iteration problem of the type un+1 = F (un), where F is the nonlinear reflection relation. We say that the PDE system is chaotic if the mapping F is chaotic as an interval map. Algebraic, asymptotic and numerical techniques are developed to tackle the cubic nonlinearities. We then define a rotation number, following J.P. Keener [11], and obtain denseness of orbits and periodic points by either directly constructing a shift sequence or by applying results of M.I. Malkin [17] to determine the chaotic regime of α for the nonlinear reflection relation F , thereby rigorously proving chaos. Nonchaotic cases for other values of α are also classified. Such cases correspond to limit cycles in nonlinear second order ODEs. Numerical simulations of chaotic and nonchaotic vibrations are illustrated by computer graphics. Received by the editors July 20, 1995 and, in revised form, October 16, 1996. 1991 Mathematics Subject Classification. Primary 35L05, 35L70, 58F39, 70L05. The first and third authors’ work was supported in part by NSF Grant DMS 9404380, Texas ARP Grant 010366-046, and Texas A&M University Interdisciplinary Research Initiative IRI 96-39. Work completed while the first author was on sabbatical leave at the Institute of Applied Mathematics, National Tsing Hua University, Hsinchu 30043, Taiwan, R.O.C. The second author’s work was supported in part by Grant NSC 83-0208-M-007-003 from the National Council of Science of the Republic of China. c ©1998 American Mathematical Society